海上OBC三维地震多波多分量采集并行模拟及应用研究

【摘要】全球陆地资源经过长期的开采,现很少有较大油气的发现。近年来海上油气储量的发现成为新增能源的重要来源。随着我国海域油气勘探程度的不断提高,勘探的难度也不断加大,还面临着一系列重大的技术挑战。海底OBC采集技术、宽方位／多方位采集以及多震源激发技术等多种新型技术被用来解决勘探中遇到的难题。目前国外海上勘探技术依然领先于国内,先进的采集技术很多已经投入到实际生产中。而针对我国海域的实际勘探情况,如何将这些先进的采集技术应用到实际生产中成为了现阶段需要解决的问题。海上外业试验成本高、耗时费力,而且得出的参数不具有完备性。而基于三维波动方程的并行数值模拟能够模拟复杂介质中波的传播特征,对各种不同的观测系统采集以及处理技术进行理论研究,为实际海上OBC采集提供理论基础和优化的采集参数。本论文针对目前海底电缆OBC采集中所采用的先进技术,主要开展三维正演模拟研究、海底OBC波场特征分析以及多震源激发与混合波场分离的研究,重点解决以下三个问题：1)针对海上气枪震源子波频带宽、主频高,海底OBC采集采用宽方位长排列的Patch观测系统而导致正演模型较大等问题,实现适用于海底OBC采集的高精度和高效率的三维正演并行算法；2)针对目前对海底OBC地震波场(特别是海底斯通利面波)传播认识存在的不足,通过正演模拟了解海底OBC采集得到的波型以及传播特征；3)实现多震源随机激发下的二维和三维海底OBC混合炮集的分离以及实际海上混合数据的分离；通过研究,本论文取得的主要成果有：1)基于波动方程有限差分数值模拟,能够模拟海底OBC采集的多波多分量数据,其中包括水听器分量、Vx速度分量、Vy速度分量以及Vz速度分量,并且基于旋度和散度求取纵波分量和转换波分量,其中分离后的纵波分量与水听器分量相似。针对海上气枪主频高、海底OBC采集通常采用宽方位长排列的Patch观测系统导致三维正演模型较大等难点,将FCT技术应用到有限差分中,正演时可以适当的增大网格大小,而FCT有限差分能够压制由网格较大引起的数值频散,取得与高阶有限差分相同的精度。2)基于MPI消息传递的并行模式可以实现三维波动方程的并行化计算,一方面将所需内存分布到不同的计算节点上,另一方面各个进程同时进行计算,提高计算效率；在并行计算过程中,可以任意选择参与并行计算的方向以及并行计算的进程个数。随着计算进程的增加,每个进程分配的计算网格点数减少,单个进程占用内存减少,正演计算时间明显减少,同时加速比随着进程数的增加而增加。但是由于串行部分的存在以及交换进程数的增加带来通信时间上的增加,故加速比增加速度逐渐变小；单个方向上并行,满足每个进程的交换进程最少(最多为2个),三个方向上并行,满足每个进程的交换信息数最少。3)海底OBC采集时,海底电缆会接收到斯通利面波。斯通利波的传播速度不仅小于海水和海底的纵波速度,同时也小于海底的横波速度,斯通利波只沿着海底界面滑行,其振幅随着与海底距离的增大而快速衰减。当海底为均匀介质时斯通利波没有频散现象,而当海底为层状介质时,斯通利波会产生频散现象。4)基于波场快照以及炮集记录波场特征分析可知,水中激发海底OBC接收的地震波场特征复杂,地震波类型很多。当震源距离海底较近时,地震波通过海底界面时主要存在的波型有海水中的直达P波、海底界面产生的反射P波、海底滑行折射P波、透射P波、透射SV波、转换折射波(PS首波)、斯通利波等。震源离海底界面较远时,斯通利波同样存在,但是能量相对来说较弱。5)多震源激发的波场是非相干波场,主要受多震源参数(震源个数、震源位置以及随机激发时间)的影响；单炮之间的距离较大时,各个炮激发的波场能明显的分开；单炮距离较小时,只要随机激发时间恰当,各个震源的波场仍然没有明显的干涉,为非相干波场；随着混合炮中震源个数的增加,地震波场也变得更加复杂。多震源混合激发得到的混合炮,经过最小二乘求逆可以得到伪分离记录,伪分离记录从共炮点道集分选到其他道集时,次震源混叠信号就会转化为随机分布的噪声。只要采用去噪的方法将随机混叠噪音去除,即可实现多震源混合炮集的分离。基于迭代的多级中值滤波去噪法,较大窗口能够去除混叠噪音,但丢失了较多的细节信息,逐渐减小中值滤波窗口大小可以提取更多的细节部分,但同时小窗口的去噪能力也在下降。本论文的创新之处主要表现在以下三个方面：1)实现了从空间X、Y和Z三个方向同时并行的三维任意偶阶空间精度的正演模拟并行算法,可以任意选择空间并行计算方向以及并行计算的进程个数；在并行算法实现过程中,引入了进程坐标系的概念,使得每一个参与计算的进程对应一个进程坐标,通过进程坐标来判断边界条件的使用；同时引入虚拟进程方便模型边界处的数据交换。2)通过数值模拟丰富了对海底OBC波场类型以及传播特征的认识,尤其是海底斯通利面波的存在以及传播特征；详细推导了海底斯通利波的表达式,并根据固液边界条件,给出斯通利波频散方程,从理论上证明了斯通利波的存在；利用Radon变换求取斯通利波的相速度频散曲线,可以得到斯通利波的传播速度。3)实现了多级中值滤波迭代去噪法进行多种二维观测系统的多震源激发的OBC混合波场分离,并在此基础上进行了三维多震源激发下的海底OBC混合炮的分离。最后将该方法应用到渤海中部石臼坨凸起区带采集到的混合数据中,经过分离后的单炮记录清晰,具有实际应用价值。
【作者】成景旺；
【导师】顾汉明；
【作者基本信息】中国地质大学，地球探测与信息技术，2014，博士
【关键词】海底电缆；多波多分量；三维并行；斯通利波；多震源激发；

【参考文献】
[1]黄海丽.基于Qt的心音分析系统的设计[D].济南大学,计算机技术,2012,硕士.
[2]邓明光.重大事件下中国股市跳跃行为特征分析[D].南京财经大学,金融学,2012,硕士.
[3]王兆翔.关于固体电解质钽电容器储存稳定性(下)[J].电子质量.2000(04)
[4]邓自立,刘玉梅.随机控制系统稳态Kalman滤波器新算法[J].自动化学报,2000,01:78-82.
[5]封雪.高中英语批判性阅读实验研究[D].河北师范大学,学科教学,2014,硕士.
[6]赖建云.论WTO环保例外权[D].湘潭大学,国际关系,2003,硕士.
[7]董珍.三次采油用Gemini表面活性剂的设计合成及性能测试[D].西北大学,无机化学,2014,硕士.
[8]李明芳.基于广义计算和分类回归树的电信业异动客户识别[D].浙江大学,2003.
[9]乔田田.马克思的文化思想及其对当代中国文化建设的启示[D].曲阜师范大学,马克思主义基本原理,2013,硕士.
[10]李欣彤.新型多彩的荧光碳量子点的制备以及生物成像研究[D].辽宁师范大学,2014.
[11]王宏斌.消息中间件技术的研究和实践[D].浙江大学,计算机应用,2004,硕士.
[12]张福旺.Al-Ti-C细化剂对A356铝合金车轮力学性能的影响[D].燕山大学,材料工程,2014,硕士.
[13]胡浩.基于时变参数Copula的ΔCoVaR度量技术[D].浙江工商大学,金融学,2012,硕士.
[14]李国樑.生物样品中药物浓度的HPLC-MS/MS测定及其药代动力学研究[D].南方医科大学,药物化学,2012,硕士.
[15]刘衡.海量多媒体数据的地理信息标注技术及其应用[D].中国科学技术大学,2014.
[16]贾海军.海华公司发展策略设计方案[D].清华大学,2004.
[17]张羽.我国上市公司高管薪酬粘性研究[D].吉林大学,企业管理,2013,硕士.
[18]徐生芹.中文逗号分类方法研究[D].苏州大学,计算机应用技术,2013,硕士.
[19]唐娅.新时期以来中国电影中的重庆元素研究[D].西南大学,传播学,2014,硕士.
[20]朱海艳.FZ医院融资租赁项目的风险管理探析[D].南京大学,工商管理（专业学位）,2013,硕士.
[21]刁朋娣.城市轨道交通网络抗毁性分析与仿真研究[D].北京交通大学,2014.
[22]兰海洋,金浩,刘维宁,孙晓静.双弹性垫板刚度对扣件减振性能影响研究[J].铁道标准设计,2015,02:9-12.
[23]烟建华,任廷革,王彤,刘晓峰,张帆.基于数据库的中医理论文献的信息化研究[J].中国中医基础医学杂志,2004,10:59-63.
[24]应晓蕊.工业以太网的实时性研究及系统设计[D].浙江大学,控制理论与控制工程,2004,硕士.
[25]李茂.中国实际负利率与高储蓄并存的原因与影响研究[D].首都经济贸易大学,国民经济学,2013,硕士.
[26]耿自才.微波氢等离子体发射光谱诊断[D].大连理工大学,应用化学,2004,硕士.
[27]马思琪.信息加工心理学取向的微博可信度研究[D].上海交通大学,新闻传播,2014,硕士.
[28]吴思廉.中学教师的反思、应付方式、自动思维及其关系的研究[D].福建师范大学,教育,2003,硕士.
[29]左前明.大采高综采工作面煤尘扩散规律及防治技术研究[D].中国矿业大学（北京）,安全技术及工程,2014,博士.
[30]王传芳.地震危险性分析中地震动衰减的不确定性校正[D].哈尔滨工业大学,结构工程,2013,硕士.
[31]李荟.中亚五国消除贫困策略研究[D].华东师范大学,外交学,2013,硕士.
[32]王英.港口物流中的流程知识挖掘研究和智能优化设计[D].北京交通大学,2014.
[33]刘环.张謇南通社会自治模式的理性思考[D].西北大学,行政管理,2014,硕士.
[34]杨伟军.基于地震预警信息的火车紧急处置系统设计与研究[D].成都理工大学,信号与信息处理,2013,硕士.
[35]刘磊.大壁虎后肢中央模式发生器（CPG）神经网络的初步研究[D].南京航空航天大学,机械设计及理论,2013,硕士.
[36]杨爱萍.加工中心人机系统分析设计及改善[D].西安建筑科技大学,机械设计及理论,2004,硕士.
[37]杨洁.当代中国社会道德的网络话语呈现[D].兰州大学,社会学,2013,硕士.
[38].2014年度国家技术发明奖获奖项目目录(通用项目)[J].中国科技奖励,2015,02:27-33.
[39]张宁,贾自艳,史忠植.使用KNN算法的文本分类[J].计算机工程,2005,08:171-172+185.
[40]刘忠义,刘崇茹,李庚银.提高直驱永磁风机低电压穿越能力的功率协调控制方法[J].电力系统自动化,2015,03:23-29.
[41]冯爱玲.GM系列直流断路器级差配合的研究与分析[D].华北电力大学（河北）,电气工程,2004,硕士.
[42]刘畅.R企业员工乐观程度实证研究[D].首都经济贸易大学,企业管理,2013,硕士.
[43]涂占新.数据挖掘方法及其应用展望[J].中南财经政法大学学报,2003,02:117-120.
[44]庞鸿飞.伦敦奥运会中美女篮进攻技术对比研究[D].北京体育大学,体育教育训练学,2013,硕士.
[45]艾龙飞.楚文化与荆楚武术发展战略研究[D].华中师范大学,民族传统体育学,2013,硕士.
[46]王柄钱.农村新型社会养老保险制度研究[D].浙江海洋学院,农业推广（专业学位）,2013,硕士.
[47]彭云飞.城市规划馆空间设计研究[D].吉林大学,艺术设计（专业学位）,2013,硕士.
[48]彭晓辉.杭州西湖水域原生动物物种多样性及缘毛目一新属/种的形态学研究[D].杭州师范大学,动物学,2007,硕士.
[49]李朝阳.马钱子在治疗风湿病中的配伍减毒与增效研究[D].河南中医学院,中医骨伤科学,2012,硕士.
[50]束文举,何林霞,蒋双双,刘双玲,郭进,陆健,高献礼.甜油滋味物质及滋味特征研究[J].食品工业科技.